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ABSTRACT 

\^ , Deep spectropolarimetric observations, obtained with the Very Large Telescope (VLT), 

are presented for two powerful radio galaxies, 0850-206 (z=1.3373) and 1303+091 
(z=1.4093). These observations cover the rest-frame wavelength range ~ 1450 -3750 A. 
New radio observations and continuum images of the same sources are also presented. 
I These galaxies are the first two observed from a complete sample of nine radio sources 

. with redshifts in the range 1.3 < z < 1.5 (selected from the equatorial sample of pow- 

2^ ' erful radio sources of Best, Rottgering & Lehnert), as part of a project aimed to 

investigate the multi-component nature of the UV continuum in radio galaxies and, 
in particular, any variations of the continuum properties with the radio source age. 
' The larger radio source of the two, 0850-206, presents a high continuum frac- 

, tional polarization, averaging 17 per cent across the observed wavelength range and 

reaching 24 per cent at rest-frame wavelengths of < 2000 A. The smaller radio source, 
_ ^ 1303-f091, shows a lower continuum polarization, averaging 8 per cent and rising to 

11 per cent for rest-frame wavelengths ^ 3000 A. For both galaxies, the position angle 
Q ' of the electric vector is generally constant with wavelength and within ^ 15° of per- 

^ I pendicular to the radio axis. Both their total flux spectra and polarized flux spectra 

C/3 i reveal the 2200 A dust feature, and comparison with dust scattering models suggests 

' that the composition of the dust in these galaxies is similar to that of Galactic dust. In 

0850-206, scattered quasar radiation dominates the UV continuum emission, with the 
, , nebular continuum accounting for no more than 22 per cent and no requirement for 

' any additional emission component such as emission from young stars. By contrast, in 

^ , 1303-f091, unpolarized radiation could be a major constituent of the UV continuum 

emission, with starlight accounting for up to ^ 50 per cent and the nebular continuum 
accounting for ^11 per cent. 

The emission-line properties of the galaxies are also studied from their total inten- 
sity spectra. Comparison of the measured emission-line ratios with both shock- and 
photo-ionization models shows that the nuclear and extended gas in these galaxies is 
mainly photoionized by the central active nucleus. 

Key words: galaxies: active — galaxies: ISM — polarization — scattering — galax- 
ies: evolution — radio continuum: galaxies — galaxies: emission lines. 



1 INTRODUCTION 

The galaxies associated with powerful extragalactic ra- 
dio sources are uniquely important for understanding the 
physics of active galactic nuclei (AGN) and for studying the 
relationship between radio source activity and the proper- 
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ties of the host galaxy and its environment. Powerful dis- 
tant radio galaxies are thought to be the progenitors of 
present day, giant ellipticals (e.g. Best, Longair & Rottgering 
1998; McLure & Dunlop 2000). However, compared with 
normal elliptical galaxies, powerful radio galaxies at z ^ 0.6 
show enhanced optical/UV continuum and line emission, 
which is generally elongated and aligned along the radio 
axis (e.g. McCarthy et al. 1987). 
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The emission-line structures usually extend large dis- 
tances from the nucleus (5 ^ 100 kpc; e.g. Tadhunter 1986; 
Baum et al. 1988); the properties of these extended emission- 
line regions (EELR) can bo greatly influenced by shocks 
resulting from interactions between the radio source struc- 
tures and the interstellar/intergalactic medium (ISM/IGM) 
(e.g. Clark et al. 1998; Solorzano-Inarrea, Tadhunter & Axon 
2001; Solorzano-Inarrea & Tadhunter 2003). Interestingly, 
the emission-line properties of a sample of z~l radio galax- 
ies have been found to evolve strongly as the radio source 
passes through the host galaxy (Best, Rottgering & Lon- 
gair 2000): large radio sources (^150 kpc) show quiescent 
'rotation-dominated' velocity profiles and ionization states 
in agreement with AGN-illumiiiatioii, while smaller radio 
sources present highly distorted kinematic profiles and over- 
all ionizations consistent with being shock-dominated. 

As for the optical/UV continuum excess observed in 
radio galaxies, several different mechanisms are known to 
contribute to it (e.g. Tadhunter et al. 2002), but the relative 
extent of their contributions remains uncertain. A popular 
early interpretation for this UV excess was recent star format 
tion, induced by the passage of the radio jet through the ISM 
of the host galaxy or by merger events linked to the radio 
source triggering (e.g. Rees 1989). Jet-induced star forma- 
tion has been shown to be feasible in numerical simulations 
(e.g. Mellema, Kurk & Rottgering 2002). Also, stellar ab- 
sorption features characteristic of OB stars have been ob- 
served in the radio galaxy 4C41.17 at z=3.8 (Dey et al. 1997). 

The detection of polarized continuum and polarized 
broad permitted emission lines in some radio galaxies (e.g. 
Antonucci & Miller 1985; di Sorogo Alighieri et al. 1989; 
Cimatti et al. 1993) suggested that light emitted anisotrop- 
ically by a hidden quasar nucleus and scattered towards us 
by dust or electrons in the ISM of the radio galaxy (Tad- 
hunter, Fosbury & di Serego Alighieri 1988; Fabian 1989) 
also makes an important contribution to the UV light. This 
is favoured by orientation-based unified schemes for radio 
sources (Barthel 1989), according to which radio galaxies 
and quasars are drawn from the same parent population 
but viewed from different angles to the line of sight, with 
the AGN in radio galaxies obscured by a surrounding dusty 
torus. Spectropolarimetric studies of a small number of radio 
galaxies at z~l have shown that the polarized emission is 
spatially extended, with the electric vector oriented perpen- 
dicular to the UV continuum emission at all wavelengths, 
and that, while the permitted Mgll 2800 omission lino is ob- 
served to be broad in polarized light, none of the narrow 
lines shows significant polarization (e.g. Cimatti et al. 1996, 
1997). These results clearly show that the origin of the polar- 
ization is scattered quasar light. More recent spectropolari- 
metric studies of a sample of radio galaxies at z~2.5 (Vernet 
et al. 2001) have shed light on the nature of the scattering 
material. These studies show indications of a continuum up- 
turn beyond 2200 A in the composite spectrum of the galax- 
ies in the sample, which is interpreted as a possible detection 
of the 2200 A dust feature, indicating that the scattering 
medium is dust. However, caution must be taken because, 
due to the high redshift of their sample, those wavelengths 
are redshifted into a region of strong sky lines. 

Together with the young stellar and scattered AGN 
components, other processes known to contribute to the UV 
continuum excess in some sources are nebular continuum 



emitted by the extended emission-line gas (Dickson et al. 
1995) and direct AGN light (Shaw et al. 1995). 

As with the emission-line gas, some properties of the 
continuum emission in radio galaxies are observed to vary 
with radio source size: smaller sources show strings of bright 
knots aligned along the radio axis, and larger sources present 
more compact nebulae with fewer bright components (Best, 
Longair & Rottgering 1996). Do the contributions of the 
different components to the UV excess also evolve over the 
radio source lifetime? For instance, a large variation is seen 
in the strengths, and to a lesser extent colours, of the po- 
larized emission in radio galaxies. Is this linked to the radio 
source size? Such questions cannot be answered by existing 
spectropolarimetric data, because targets were generally se- 
lected on the basis of previously detected polarization in 
imaging observations, an interesting UV morphology or an 
ultra-steep radio spectrum, making the samples studied bi- 
ased or incomplete. 

To remedy this and address many of the issues raised 
above, such as the nature of the scatterers at high redshift 
and the variation with radio source size of the different con- 
tributions to the UV excess, we have begun a programme to 
make deep spectropolarimetric observations of a complete 
sample of nine radio galaxies with redshifts in the range 
1.3 <z< 1.5 and a wide range in radio sizes. In this paper 
we present the observations of the first two galaxies from 
the sample. 

The paper is organized as follows. In Section 2, the selec- 
tion criteria of our sample are described. Section 3 contains 
the details of the observations, data reduction and analysis. 
The results are presented in Sections 4 and 5. Discussion of 
the results follows in Section 6. Summary and conclusions 
arc provided in Section 7. Throughout this paper values 
of the cosmological parameters of Ho = 65kms~^ Mpc~^, 
Q,M = 0.3 and Qa =0.7 are assumed. 



2 SAMPLE SELECTION 

The sources of our full sample comprise all radio galax- 
ies from the equatorial sample of powerful radio sources 
of Best, Rottgering & Lehnert (1999)^, with redshifts in 
the range 1.3<z<1.5. The redshift range was chosen to 
be low enough to allow good S/N, but also sufficiently high 
to ensure that the rest-frame wavelengths studied (~ 1400- 
3800 A) allowed to search for features in the dust scattering 
cross-section around 2200 A, as well as stellar absorption 
lines, such as Silll 1417, SV 1502 and NIV 1720. Further, this 
wavelength range covers the Hell 1640 emission line, which 
can be used as a replacement for the Balmer lines to derive 
the amount of nebular continuum emission. Our sample also 
includes sources with a wide range of radio sizes, from ~ 10 
to 400 kpc, in order to investigate variations of the emission- 
line and contirmum properties with the radio source age. 

The results obtained from the first two sources observed 
from our sample are presented in this paper. Properties of 
these two sources are listed in Table 1. 

t Following the additional identifications and redshifts of Best, 
Rottgering & Lehnert (2000) and Best et al. (2003), this radio 
source sample Is now 100 per cent spectroscoplcally complete. 



© 0000 RAS, MNRAS 000, 000-000 



VLT spectropolarimetry of powerful radio galaxies 3 



Source 


z 


Drad 


■P-^rad 


Eb-v 






(kpc) 


(dcg) 


(mag) 


0850-206 


1.3373±0.0018 


118 


29 


0.242 


1303+091 


1.4093±0.0012 


73 


131 


0.030 



Table 1. Properties of the two radio galaxies from our sample 

presented in this paper. The redshift of each source, derived from 
the data presented in this paper, is given in column 2. Column 
3 gives the projected extent of the radio source along the radio 
axis (Best, Rottgering & Lehncrt 1999). The position angle of the 
radio axis for each source (Best et al. 1999) is listed in column 
4. The last column gives the colour excess Eb_v for each source 
taken from NASA/IPAC Extragalactic Database (NED). 

3 OBSERVATIONS, DATA REDUCTION AND 

ANALYSIS 

3.1 VLT observations 

Optical imaging and spectropolarimetric observations were 
carried out on the 8.2-m ESQ Very Largo Telescope (VLT) 
Antu (UTl) at Paranal (Chile) on the nights of 2000 Jan- 
uary 31 and February 1. The imaging spectrograph FORSl 
was used in both direct imaging (IMG) and multi-object 
spectropolarimetry (PMOS) modes together with the Tek 
2048x2048 CCD, which provides a spatial scale of 0.2 axe- 
sec pixel" ^. Observations were made under good seeing con- 
ditions, averaging 0.5 arcsec and 0.7 arcsec for the two 
nights, respectively. The first night was clear, but part of 
the second night was non-photometric. Details of the obser- 
vations are presented in Tables 2 and 3. 

3.1.1 Optical imaging 

Imaging observations were made through the broad-band 
filters V_BESS-F35 for 0850-206, and I_BESS-f37 for both 
0850-206 and 1303+091. 

The reduction of the imaging data was performed using 
standard tasks within the IRAF software package, following 
the usual steps. The bias was first subtracted and then the 
images were flat-fielded using sky flats. Photometric calibra- 
tion was achieved by using observations of the photometric 
standard stars PC 0231+051 and PC 1323-086 taken during 
the nights. The Johnson photometric magnitudes, measured 
through 2- and 4-arcsec diameter apertures, are given in Tar 
ble 2. 

3.1.2 Optical spectropolarimetry 

The grism GRIS_300V was used for the spectropolarimet- 
ric observations, providing a spectral dispersion of approxi- 
mately 2.6 Apixel"^. For the source 0850-206 a 1. 4-arcsec 

slit was oriented along the axis defined by the position of op- 
tical objects, which is misaligned by approximately 35° with 
the radio axis [see Fig. 1 (left)]. In the case of 1303+091, a 
1-arcsec slit was oriented parallel to the radio axis [see Fig. 6 
(left)]. Note that the imaging spectrograph FORS has a com- 
pensator that corrects for atmospheric dispersion. Four sets 
of observations were made for each source (two sets on each 
night), each set consisting of four exposures with the half- 
wave plate position angle rotated to 0°, 22°5, 45° and 67°5, 



Source Filter Acen texp magnitude Airm. 

(A) (s) 2" 4" 

0850-206 V3ESS 5540 300 22.98 22.73 1.024 

0850-206 I.BESS 7680 300 21.79 21.39 1.018 

1303+091 I.BESS 7680 300 21.83 21.40 1.724 



Table 2. Log of the VLT optical imaging observations. Photo- 
metric magnitudes, measured through 2- and 4-arcsec diameter 
apertures, are also given. 



consecutively. 0850 206 was observed in three sets of four 
1350 s exposures and one set of four 900 s exposures, with 
a total integration time of 5.5 hr. 1303+091 was observed 
in two sets of four 1350 s, one set of four 900 s and one set 
of three 1150s (the 67°5 exposure was lost due to instru- 
ment problems), making a total integration time of nearly 
5hr. The lEist two sets of 1303+091 were carried out in non- 
photometric conditions. 

The reduction of the spectropolarimetric data was per- 
formed using IRAF in the following way. The frames were 
first bias-subtracted and then fiat-fielded using lamp flats 
taken through the same optics used in the source and star 
frames. After removing cosmic-ray events, the different ex- 
posures for the same half-wave plate positions were com- 
bined for each source. In the case of 1303+091, appro- 
priate weights and scales were applied to account for the 
non-photometric conditions in which this source was ob- 
served during the second night. One-dimensional spectra, 
with identical aperture widths for the ordinary (o) and ex- 
traordinary (e) rays for each of the four half-wave plate ori- 
entations, were then extracted. The wavelength calibration 
was performed using He, HgCd and Ar arc lamps. After 
that, o- and e-rays were resampled to the same linear spec- 
tral dispersion, in order to calculate the polarization over 
identical spectral bins. The data were then binned, and all 
eight one-dimensional spectra were combined to determine 
the Stokes parameters I, Q and U for each galaxy, using 
the procedures described in Vernet (2001; see also Vernet 
et al. 2001), which are based on the method outlined in Co- 
hen et al. (1995). Observations of polarized (HD245310 and 
HD111579) and unpolarized (HD64299 and HD94851) stan- 
dard stars (Turnshek et al. 1990; Schmidt, Elston & Lupie 
1992) were taken in order to check and calibrate the instru- 
mental polarization and the position angle zero-point offset. 
Flux calibration was provided by observations of the spec- 
trophotometric standard stars Feige24 and G163-50. The 
atmospheric extinction correction was achieved by using the 
extinction curve of La Silla Observatory, since no suitable 
data for Paranal arc available. 

No second-order filter was used in these observations, 
but instead the polarized standard star HD245310 was ob- 
served both with and without the GG435 filter, which blocks 
the light below 4200 A, to estimate the effects of any second- 
order blue light. Residual second-order light, which only af- 
fects the spectral region with A > 7000 A, was found to be 
very low. In addition, the continuum polarization of the star 
HD245310 was derived for both sets of observations (with 
and without filter), and it was found that it did not vary 
in any significant way. Thus, any contamination by second- 
order blue light does not affect the conclusions of this paper. 
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Source 


texp 


Slit PA 


Grism 


A Range 


Pixel Scale 


Resolution 


Slit Width 


Airmass 




(s) 


(deg.) 




(A) 


(arcscc/A) 


(A) 


(arcscc) 




0850-206 


19800 


175 


GRIS.300V 


3460 - 8792 


0.2/2.607 


14.66±0.10 


1.4 


1.19 


1303+091 


17850 


131 


GRIS.300V 


3281 - 8583 


0.2/2.593 


11.27±0.07 


1.0 


1.36 



Table 3. Log of the VLT optical spectropolarimetric observations. The spectral resolution given is that of the observed frame. 



Source 


Freq. 


Array 


Obs. Date 


texp 




(MHz) 


(Config.) 


(d/m/y) 


(min) 


0850-206 


8460 


A 


09/12/2000 


56 






B 


20/05/2001 


16 


0850-206 


4710 


A 


09/12/2000 


38 






B 


20/05/2001 


15 


1303+091 


8460 


A 


11/12/2001 


55 






B 


20/05/2001 


16 


1303+091 


4710 


A 


11/12/2001 


33 






B 


20/05/2001 


15 



Table 4. Log of the VLA radio observations. 

For the analysis of the total flux spectra, the two- 
dimensional frames of the o- and e-rays were first cor- 
rected for any distortion along the spectral direction, us- 
ing the Starlink FIGARO package, and then summed. Af- 
ter that, one-dimensional spectra were extracted, and then 
analysed using the Starlink DIPSO spectral analysis package. 
Emission-line fluxes, velocity shifts and linewidths (FWHM) 
were obtained by Gaussian fitting of the emission-line pro- 
files. The measured linewidths were corrected for the spec- 
tral resolution of the instrument via Gaussian dcconvolution. 
The instrumental resolution was measured from night-sky 
lines and arc lines, and is listed in Table 3. 

The total intensity spectra, emission-line fiuxes and ra- 
tios have been corrected for Galactic reddening using the ex- 
tinction curve from Seaton (1979) and values of colour excess 
Eb-v from NASA/IPAC Extragalactic Database (NED; see 
Table 1). 

3.2 Radio observations 

New radio observations of all nine sources in the spectropo- 
larimetric sample were carried out using the Very Large Ar- 
ray radio interferometer. Each source was observed using the 
A-array configuration at both 8.5 and 4.7 GHz frequencies, 
mapping both the total intensity and polarization proper- 
ties. Sources larger than 10 arcsec in angular extent were 
additionally mapped usiug the B-array configuration in or- 
der to provide sensitivity to the extended structures. Full de- 
tails of all of these observations will be provided elsewhere 
(Best et al., in prep.). Here, only the total intensity data 
for 0850-206 and 1303+091 are considered, to compare the 
radio properties with the spectropolarimetric and imaging 
data and to use the radio spectral index information to lo- 
cate the radio cores. 

Details of the observations of these two sources arc pro- 
vided in Table 4. Data reduction for both sources followed 
the standard procedures of CLEANing and self-calibration 
within the AIPS package (e.g. see Best et al. 1999). 8.5 GHz 



total intensity maps were made at full angular resolution 
(« 0.3 arcsec), with rms noise levels of 22/iJy. Maps were 
also made at lower angular resolutions matching those of the 
4.7 GHz data (0.45 arcsec for 1303+091; 0.67 by 0.45 arcsec 
for 0850-206), by weighting down the outer regions of the 
UV-plane. These matched-resolution 8.5 GHz and 4.7 GHz 
maps were aligned and combined to produce a radio spec- 
tral index map of each source. 

3.3 Registration of the radio and optical images 

The optical and radio maps of the two galaxies were aligned 
by assuming that the continuum centroid of the optical 
emission was coincident with the location of the radio core, 
adopted to be the radio knot with the flattest spectral index 
[sec Figs. 1 (right) and 6 (right)]. As a check, astrometry was 
also performed on the optical images by using several unsat- 
urated stars present in the observed fields and in the APM 
(Maddox et al. 1990) and US Naval Observatory catalogues. 



4 0850-206 

4.1 The continuum emission 

Fig. 1 (left) presents the VLT image of the radio galaxy 
0850-206 in the V band, overlaid with contours of the radio 
emission at 8.5 GHz. The position of the slit for the spec- 
tropolarimetric observations is also shown in the figure. The 
central wavelength of the V-filter in the rest frame of the 
galaxy was 2370 A, with a bandwidth of 440 A (rest-firame). 

Fig. 2 shows the spatial profile of the line-free contin- 
uum emission (solid line) along the slit for 0850-206, for 
the rest-frame wavelength ranges 1700 1850 A (left) and 
2900-3050 A (right). It appears that the total continuum 
emission in the UV end is slightly weaker than that at longer 
wavelengths. 

4.1.1 Polarization of the continuum emission 

The spectropolarimetric results of 0850-206 are presented 
in Fig. 3. From top to bottom the following are plotted: the 
total flux spectrum, the percentage polarization and the po- 
sition angle of the electric vector. A 2 arcsec-wide (~ 18 kpc) 
aperture centred on the continuum centroid was used for the 
extraction of ID spectra for the o- and e-rays. 

The continuum fractional polarization of this source is 
observed to increase gradually towards shorter wavelengths, 
with an average value of 17 per cent across the observed 
wavelength range of the spectrum. The highest polarization 
is measured in the bin longward of Hell 1640, with a value 
of 24.2I4 gper cent. Note that the polarization measured in 
the bin shortward of [OH]3727, with a value of 6 per cent, 
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Figure 1. Left: V-band image of 0850—206 overlaid with contours of the radio emission at 8.5 GHz. Tlie values of the radio eontours are 
1x10"** Jybeam-l X (-1, 1, 2, 4, 8, 16, 32, 64, 128, 256, 512). The position of the slit is also indicated on the left plot (the length of 
the slit is arbitrary). Right; Radio spectral index, a (greyscale, defined as Si^ oci'~"), between 8.5 and 4.7 GHz for 0850—206. Overlaid 
are contours of radio emission at 4.7GHz, with contour levels of 1.7x lO"** Jy beam"* x (-1, 1, 2, 4, 8, 16, 32, 64, 128, 256, 512). The 
bright radio knot at RA = 08h53m00.27s Dec = -20° 47' 27.6" has the flattest spectral index, aRiO.4 which is typical of radio cores at 
these redshifts (e.g. Athreya et al. 1997). 



may be affected by noise due to the strong sky lines in that 
spectral region. 



The position angle of the electric vector is generally con- 
stant with wavelength, with an average value of 103°, which 
is within ^ 15° perpendicular to the radio axis (PArad ~ 29°, 
defined by the position of the two most distant radio 
hotspots). For shorter rest-frame wavelengths than 1900 A, 
the position angle seems to decrease slightly down to a 
value of ~ 90° . Previous studies of high-z radio galaxies 
show a better perpendicularity between the electric vec- 
tor and the UV continuum axis than with the radio axis 
(e.g. di Serego Alighieri, Cimatti & Fosbury 1993; Cimatti 
et al. 1993, 1994; di Serego Ahghieri et al. 1996; Vernet et al. 
2001). There are no rest-frame UV images of 0850-206 avail- 
able to measure the UV axis with precision, but the UV im- 
age published in this paper [ Fig. 1 (left)] , although not very 
deep, suggests that the UV axis is 12° ±5°, which is perpen- 
dicular to the electric vector orientation. It should also be 
noted that previous studies of low-z narrow-line radio galax- 
ies have shown that when there is a hidden nucleus and a 
polarization fan, the measured position angle of the electric 
vector can depend on the slit orientation (Cohen et al. 1999). 



4.2 The line emission 

4-2.1 Emission-line structure 

2D spectra of some of the main emission lines observed for 
0850-206 are shown in the top panel of Fig. 4. These are 
Hell 1640, CIII]1909, CII]2326, [NeIV]2425, [NeV]3426 and 
[Oil] 3727. The spatial variation of the flux of the brightest 
emission line ([Oil] in this case) along the slit was deter- 
mined and it is presented in the same figure (bottom panel 
- top left). The spatial profile weis derived by first extract- 
ing a spatial slice perpendicular to the dispersion direction 
and centred on the [Oil] emission line, with a spectral width 
of 60 A. Then, to correct for the continuum contamination, 
a spatial slice extracted from a line-free region of the spec- 
trum, adjacent to the emission line and scaled to have the 
same spectral width as that of the emission-line slice, was 
subtracted from the first one. 

The detected [Oil] emission extends approximately 6 
arcsec (~54kpc) along the slit. The line emission distribu- 
tion is observed to be asymmetric with respect to the nu- 
cleus, extending further to the south than to the north. 

4-2.2 Emission-line kinematics 

Kinematic information of the emission-line gas in 0850 206 
was obtained by fitting single Gaussians to the profile of 
[Oil] 3727, which was the emission line with the highest 
signal-to-noise ratio in the observed spectrum. To try to 
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Figure 2. Spatial profiles of tlic linc-frcc continuum emission (solid line) along the slit for 0850—206, between 1700 — 1850 A rest-frame 
(left) and between 2900 — 3050 A rest- frame (right). The dashed line shows the nebular continuum profile for the same wavelength ranges 
(see Section 6.1). South is to the left. 
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Figure 3. Spectral and polarization properties of the radio galaxy 0850-206. Prom top to bottom: the total flux spectrum, the observed 
fractional polarization and the position angle of the electric vector. Filled circles and crosses indicate respectively continuum bins and 
emission line bins with their underlying continuum. 



match the seeing of the observations, one-dimensional spec- 
tra with an aperture of 3 pixels (0.6 arcsec) and centred in 

steps of 2 pixels were extracted from the 2D total inten- 
sity spectrum of the galaxy. In the outermost regions of the 
galaxy the aperture width of the extracted ID spectra was 
greater than 3 pixels in order to increase the S/N of these 
regions. 



The spatial variation along the slit of the velocity cen- 
troids and hnewidth of the [Oil] 3727 line for 0850-206 are 
shown in Fig. 4 (bottom panel - top right, bottom right). 

The velocity shifts, which arc referred to the velocity 
of the emission-line gas at the continuum centroid, are ob- 
served to vary smoothly along the position of the slit, with 
an overall velocity amplitude of Aw ~ 450 kms~^. The ve- 
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Figure 4. Top panel: 2D spectra of the main emission lines for 0850-206. Each spectrum covers 9.6 arcsec along the slit, and 30 A 
(rest-frame) in the spectral direction. The regions defined along the slit are indicated by the solid line at the right of the panel. Bottom 

panel: Emission-line spatial analysis along the slit for 0850-206. Top left: Variation of the [011]3727 flux. The spatial regions defined 
along the slit are indicated by a solid line at the top of the plot. Bottom left: Spatial variation in the CIIIJ1909 / CII]2326 line ratio. Top 
right; Spatial variation in the [OII]3727 velocity centroids. Bottom right: Spatial variation of the [OII]3727 linewidth. South is to the left 
in the four plots. 



locity curve is, however, highly asymmetric with respect to 
the nucleus; while the emission-line gas at the north end, 
about 1 arcsec from the centre, appears to have a velocity 
comparable to that of the nucleus, the gas at ~ 3 arcsec to 
the south of the nucleus is redshifted by about -1-400 kms~^. 

The instrumentally-corrected linewidths (FWHM) of 



the [Oil] line are found to range between ~300 and 800 
kms~^ along the location of the slit. The broadest emission 
lines, with FWHM'~ 750kms~^, are located about 1 arc- 
sec south of the nucleus. The nucleus itself presents [Oil] 
linewiths of FWHM ~ 620 kms"^ 
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Figure 5. Integrated spectra for the different emission-line re- 
gions in 0850-206 defined along the slit: northern EELR (top; 
1.2x1.4 arcsec^ aperture centred at 1.1 arcsec N of the continuum 
centroid), nucleus (middle top; 1.2x1.4 arcsec^ aperture centred 
at 0.1 arcsec S of the continuum centroid), southern EELR 'SI' 
(middle bottom; 1.2x1.4 arcsec^ aperture centred at 1.3 arcsec 
S of the continuum centroid) and southern EELR 'S2' (bottom; 
2.2x 1.4 arcsec^ aperture centred at 3.0 arcsec S of the continuum 
centroid) . 



4-2.3 Emission-line spectra 

The total intensity spectrum of 0850 206 (corresponding to 
an aperture of 2 arcsec centred on the nucleus) is shown in 
Fig. 3 (top). Also, integrated spectra of four different regions 
defined along the slit (nuclear region, southernmost 'S2', 
southern 'SI' and northern EELR) arc presented in Fig. 5 
(sec the caption of the figure for details of the apertures). 
The emission-line fluxes normalized to the CII] 2326 flux for 
each region are listed in Table 5. CII] 2326 was chosen for 
the normalization of the line fluxes because it is a common 
line in the observed spectra of both 0850 206 and 1303+091, 
whose emission is also detectable in the extended regions of 
the galaxies. 

The spatial variation of the CIII]1909 / CII]2326 Une rar 
tio along the location of the slit for 0850-206 is shown in 

Fig. 4 (bottom panel - bottom left). The emission- line gas 
with the highest ionization state is found to be at ~ 0.5 arc- 
sec north of the nucleus, with a line ratio CIII] / CII] ~ 3.6. 
The nucleus itself presents a moderate ionization state in 
comparison with the rest of the regions observed. The gas 
showing the lowest ionization, with CIII] / CII] ~ 1.3, is lo- 
cated at ~ 2 arcsec south of the nucleus, which coincides with 
a region of relatively broad emission lines (FWHMjou] ~ 660 
kms"^). 

4-2.4 Polarization of the emission lines 

In view of the unified schemes for radio sources, it is ex- 
pected that the forbidden lines are unpolarized, since they 
are emitted isotropically in the narrow-line region outside 
the obscuring torus (unless resonantly scattered). On the 
other hand, the broad permitted lines, emitted in the broad- 
line region inside the torus, are expected to be polarized. 

Fig. 3 (middle plot) shows the observed fractional po- 
larization for coritirmum and cmission-linc bins in 0850 206, 
for an spatial aperture of 2 arcsec centred on the contin- 
uum centroid. It can be seen that the polarization of the 
emission lines (with the underlying continuum) is lower than 
the adjacent pure continuum bins. After subtracting the un- 
derlying contirmum, the polarization of the emission lines 
with high enough S/N was measured, obtaining the fol- 
lowing: Pciii] =5.114 5 per cent, P[Neiv] =3.7±2.5 per cent, 
P[NeV] = 4.1±2.9 per cent and P[oii] = 3.6±1.1 per cent. The 
position angles of the electric vector for these lines range 
from 112° to 135°. These low polarizations (if real) could be 
due to transmission through magnetically aligned dust in ei- 
ther the radio galaxy itself or the ISM of our Galaxy: using 
the Galactic relation Pism<9.0Eb-v (Scrkowski, Mathow- 
son & Ford 1975), the interstellar polarization due to our 
Galaxy could be up to ~2.2 per cent. However, the fact 
that the electric vector position angles of these lines are on 
average almost perpendicular to the radio axis of 0850-206 
(PArad ~ 29°), suggests that the low polarization observed in 
these omission lines could be duo to scattering or residuals 
of the continuum subtraction. 

The continuum-subtracted polarization of the permit- 
ted Mgll fine was also measured, obtaining PMgii = 4.6^4 ;[! 
per cent and an electric vector position angle of 148° ±24°. 
The errors in the polarization measurements of this line 
are bigger and thus no significant conclusions can be drawn 
about its polarization properties. 
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0850-206 S2 EELR SI EELR Nucleus N EELR 

Line (3".0 S) (1".3 S) (0".l S) (0".9 N) 



HeII1640 <3.28 <1.06 2.56±0.21 5.22±0.95 

NIII] 1749 — — 0.32±0.13 1.07±0.46 

Silll] 1882 — — 0.18±0.06 — 

Silll] 1892 — — 0.25±0.06 — 

CIII] 1909 1.17±0.43 1.37±0.26 2.46±0.17 2.92±0.52 

CII]2326 1.00±0.23 1.00±0.13 1.00±0.06 1.00±0.16 

[NeIVl2425 1.21±0.46 1.50±0.28 2.95±0.18 2.19±0.39 

[Oil] 2472 — — 0.23±0.05 — 

[MgVII]2510 — — 0.33±0.06 — 

[MgVII]2629 — — 0.21±0.06 — 

HcII2733 — — 0.53±0.07 — 

[MgV]2783 — — 0.28±0.05 — 

MgII2800 — — 1.92±0.14 2.53±0.44 

HcI2829 — — 0.14±0.04 — 

OIII3133 — — 0.42±0.06 0.64±0.17 

Hell 3203 — 0.55±0.16 0.52±0.06 — 

[NeV]3346 — 1.20±0.24 1.80±0.13 1.12±0.28 

[NeV]3426 <1.08 3.69±0.53 4.28±0.26 1.96±0.41 

[Oil] 3727 10.88±2.56 18.43±2.44 11.50±0.69 9.54±1.55 



Oil] 2326 flux 

(10-1^ ergs-icm-2) 2.42±0.55 3.11±0.41 7.77±0.45 2.25±0.35 



Table 5. Emission-line integrated fluxes (corrected for Galactic reddening and normalized to the Cll]2326 flux) for the different regions 
in 0850-206. The details of the apertures are in the caption of Fig. 5. The quoted errors correspond to the line fitting errors. 



5 1303+091 

5.1 The continuum emission 

Fig. 6 (left) shows the VLT I-band image of the radio galaxy 
1303+091, overlaid with contours of the radio emission at 8.5 
GHz. The position of the slit for the spectropolarimetric ob- 
servations is indicated in the figure. The central wavelength 
of the I-filter in the rest frame of the galaxy was 3188 A, 
with a bandwidth of 573 A (rest-frame) . 

Fig. 7 shows the spatial profile of the line-free contin- 
uum emission (solid line) along the slit for 1303-1-091, for 
the rest-frame wavelength ranges 1700 - 1850 A (left) and 
2900 - 3050 A (right). It can be seen that, unlike 0850-206, 
1303-1-091 presents a stronger continuum emission at the UV 
end than at longer wavelengths. 



5.1.1 Polarization of the continuum emission 

Fig. 8 presents the spectropolarimetric results of 1303-1-091, 

showing from top to bottom: the total flux spectrum, the 
percentage polarization and the position angle of the elec- 
tric vector. The aperture used for the extraction of the ID 
spectra for the o- and e-rays was 2 arcsec (~ 18 kpc) wide 
and centred on the continuum centroid. 

Compared with 0850-206, 1303+091 presents an over- 
all lower continuum fractional polarization, with an average 
value of 8 per cent across the observed wavelength range of 
the spectrum. The highest polarization is found in the bin 
longward of [NeV]3426, with a value of 11.5lg | per cent. 
The polarization then, by contrast with 0850 206, is ob- 
served to decrease slightly towards shorter wavelengths. A 
dip is found in the bin of ~ 2000 A (rest frame) with a value 
of 4.8±1.0 per cent. For shorter wavelengths than ~ 1900A, 



the polarization is observed to increase again up to a value 
of 7.5^2 5 per cent. 

The position angle of the electric vector is fairly con- 
stant with wavelength, with an average value of 25° , which 
is within ~ 15° perpendicular to the radio axis, defined 
by the position of the two most distant radio hotspots 
(PArad ~ 131°), and within 7° perpendicular to the axis 
defined by the radio core and the south-eastern inner ra- 
dio knot (PAknot ~ 122°). For longer rest frame wavelengths 
than --^3300 A, the position angle of the electric vector has 
values of 64°±10° and -5+14° for the two separate bins con- 
sidered; however these values are not likely to be significant 
since that region coincides with the presence of strong sky 
emission lines, and thus it is very noisy. No comparison could 
be made between the electric vector orientation and the UV 
continuum axis because there are no rest-frame UV images 
available of 1303+091 to derive the UV axis. 



5.2 The Hne emission 

5.2.1 Emission-line structure 

2D spectra of some of the main emission lines observed for 
1303+091 are shown in the top panel of Fig. 9. These are 
CIV 1548, Hell 1640, CIII]1909, CII]2326, [NeIV]2425 and 
[NeV]3426. The spatial profile along the radio axis of the 
emission line with the highest signal-to-noise ratio in both 
nuclear and extended regions (CIII] in this case) was ob- 
tained, and it is presented in the same figure (bottom panel 
- top left). This profile was obtained in the same way as 
described for 0850-206 in Section 4.2.1, but for 1303+091 
the spectral width of the CIII] slice weis 40 A. 

The detected CIII] emission extends about 3.5 arcsec 
(~ 32 kpc) along the radio axis. The line emission is asym- 
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Figure 6. Left: I-band image of 1303+091 overlaid with contours of the radio emission at 8.5 GHz. The values of the radio contours 
are 7x 10"^ Jybeam-l X (-1, 1, 2, 4, 8, 16, 32, 64, 128, 256, 512). The position of the slit is also indicated (the length of the slit is 
arbitrary). Right: Radio spectral index, a (greyscale, defined as oc between 8.5 and 4.7 GHz for 1303+091. Overlaid are contours 

of radio emission at 4.7GHz, with contour levels of 1.5 x 10"* Jybeam"! x (-1, 1, 2, 4, 8, 16, 32, 64, 128, 256, 512). The bright radio 
knot at RA = 13h 05m 35.93s Dec = +08° 55' 14.9" has the flattest spectral index, a 0.5 which is typical of radio cores at these redshifts 
(e.g. Athreya et al. 1997). 



metrically distributed with respect to the nucleus; in addi- 
tion to the central nuclear peak, there is a fainter peat at 
about 1 arcsec (~ 9 kpc) to the north-west of the nucleus, 
close to the radio knot in the north lobe where the radio 
source appears to bend. 



5.2.2 Emission-line kinematics 

Kinematic information for the emission-line gas in 1303+091 
was obtained following the same procedure as for 0850-206 
(see Section 4.2.2). In this case, the emission line used for 
the analysis was the CIII] 1909 Une, which has the highest 
S /N in the observed spectrum and a similar ionization state 
to [Oil] 3727 used for 0850-206. 

The spatial variation of the velocity centroids and 
linewidth of CIII] along the radio axis of 1303+091 are pre- 
sented in Fig. 9 (bottom panel - top right, bottom right). 

The velocity shifts arc referred to the velocity of the 
emission-line gas at the continuum centroid. By contrast 
with 0850-206, the overall velocity amplitude around the nu- 
cleus along the radio axis of 1303+091 changes by At; 100 
kms^^ over ~3 arcsec. 

The CIII] instrumentally-corroctod linewidths 
(FWHM) along the radio axis of 1303+091 are observed 
to vary between ~600 kms~^ at 2 arcsec north-west of the 
nucleus and 900 km s~^ at 1 arcsec south-east of the nucleus. 
The CIII] linewidth at the location of the continuum cen- 
troid is found to be FWHM ~ 800 kms~^. It is interesting 
that the highest linewidths are measured at the location of 
the south-eastern inner radio knot, and that the narrowest 
lines arc found outside the radio structure of the northern 
radio lobe (see Fig. 6). This is expected if the broadening 
of the emission lines is due to interactions with the radio 
structures. 



5.2.3 Emission-line spectra 

The total intensity spectrum of 1303+091 (corresponding 
to an aperture of 2 arcsec centred on the nucleus) is shown 
in Fig. 8 (top). In addition, integrated spectra of three re- 
gions defined along the radio axis (nuclear region, south- 
eastern and north-western EELR) arc presented in Fig. 10 
(see the caption of the figure for details of the apertures). 
The emission-line fluxes normalized to the CII] 2326 flux (see 
Section 4.2.3) for each region are listed in Table 6. 

The spatial variation of the CIII] 1909 / CII] 2326 line ra- 
tio along the radio axis of 1303+091 is shown in Fig. 9 (bot- 
tom panel - bottom left). The nucleus presents the high- 
est ionization state along the radio axis, with a line ratio 
CIII] / CII] ~ 2.5. There may be a second peak in the ion- 
ization state at ~ 1.2 arcsec north-west of the contirmum 
centroid. The emission-line gas with the lowest ionization is 
located at 1 arcsec south-east of the nucleus, with a ra- 
tio CIII] / CII] ~ 1.4. Interestingly, this region of low ioniza- 
tion presents the broadest emission lines (FWHMcm] ~ 950 
kms^^), and coincides with the location of a radio knot (see 
Fig. 6). 



5.2.4. Polarization of the emission lines 

Fig. 8 (middle plot) shows the observed fractional polar- 
ization for continuum and emission-line bins in 1303+091, 

for a spatial aperture of 2 arcsec centred on the contin- 
uum centroid. It can be noticed that the polarization of the 
emission lines (with the underlying continuum) is similar 
or lower than the adjacent pure continuum bins. After sub- 
tracting the underlying continuum, the polarization of the 
stronger emission lines was measured. The permitted lines 
where found to be significantly polarized, Pciv = 9.4^4;^ per 
cent, Pneii = 11.714 3 P^'' '^^'^^ PMgii = 5.612.4 P^'' cent 
(note that these values are for the entire emission line: there 
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Figure 7. Spatial profiles of the line-free continuum emission (solid line) along the radio axis of 1303+091, between 1700 — 1850 A rest- 
frame (left) and between 2900 — 3050 A rest- frame (right). The dashed line shows the nebular continuum profile for the same wavelength 
ranges (see Section 6.1). South-east is to the left. 
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Figure 8. Spectral and polarization properties of the radio galaxy 1303-f091. Prom top to bottom: the total flux spectrum, the observed 
fractional polarization and the position angle of the electric vector. Filled circles and crosses indicate respectively continuum bins and 
emission line bins with their underlying continuum. 



was insufficient S/N to attempt to separate any broad and 
narrow components). The low Galactic reddening towards 
1303-1-091 means that the ISM in our Galaxy makes a negli- 
gible contribution to the polarization (Pism< 0.27 per cent; 
Serkowski et al. 1975). The position angles of the electric 
vector for CIV and Hell are 42° ±13° and 36° ±11°, respec- 
tively, which are almost perpendicular to the radio axis of 



1303-1-091 (PArad ~ 131°), suggesting that the polarization 
of CIV and Hell (if real) is dominated by scattering. The 
electric vector orientation for Mgll is, however, found to be 
167° ±13°. This value is far from perpendicular to the radio 
axis, and it is probably affected by the low S/N of the Mgll 
line. 

The polarization of the stronger forbidden and semi- 
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Figure 9. Top panel: 2D spectra of the main emission lines for 1303+091. Each spectrum covers 9.6 arcsec along the radio axis, 
and 29 A (rest-frame) in the spectral direction. The regions defined along the slit are indicated by the solid line at the right of the 
panel. Bottom panel: Emission-line spatial analysis along the radio axis of 1303-1-091. Top left: Spatial variation of the Clll] 1909 
flux. The spatial regions defined along the slit arc indicated by a solid line at the top of the plot. Bottom left: Spatial variation in the 
Clll] 1909 / Cll]2326 line ratio. Top right: Spatial variation in the CIII] 1909 velocity centroids. Bottom right: Spatial variation of the 
cm] 1909 linewidth. South-east is to the left in the four plots. 

forbidden lines was also measured, finding: Pciii] = 3.6^2^4 6 DISCUSSION 

per cent, P[Neiv] —'i-^-t'^, per cent and P[NeV] =6.7l|;| per 

cent, all of which are consistent with zero. 6.1 The origin of the UV continuum 

There is evidence that at least four mechanisms contribute 
to the UV excess observed in powerful radio galaxies [e.g. 
Tadhunter et al. (2002) and references therein]; these are: 
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Figure 10. Integrated spectra for the different emission-line re- 
gions in 1303+091 defined along the radio axis: north-western 
EELR (top; 1.8x1.0 arcsec^ aperture centred at 1.4 arcsec NW of 
the continuum centroid), nucleus (middle; 1.0x1.0 arcsec'^ aper- 
ture centred on the continuum centroid) and south-eastern EELR 
(bottom; 1.6x1.0 arcsec^ aperture centred at 1.3 arcsec SE of the 
continuum centroid). 



scattered AGN light, young stellar population, nebular con- 
tinuum and direct AGN light. In this section the relative 
contribution of the different components to the UV contin- 
uum of 0850-206 and 1303-1-091 is investigated. 

Direct AGN light is found to contribute significantly to 
the UV continuum of a source only when broad permitted 
lines are detected in the intensity spectrum. Such sources are 
likely to be partially obscured, or low luminosity, quasars 
and they usually present low levels of polarization of or- 
der ^ 4 per cent (e.g. Stockman, Angel & Miley 1979; An- 
tonucci 1984). In the cases of 0850-206 and 1303-^091, no 
broad components are detected in the profiles of the permit- 
ted emission lines in their spectra; also, for both galaxies the 



1303-1-091 SE EELR 

Line (1".3 SE) 

CIV 1548 2.10±0.33 

Hell 1640 1.59±0.36 

OKI] 1666 — 

[MgVI] 1806 — 

Silll] 1882 — 

Silll] 1892 — 

cm] 1909 1.08±0.13 

Nil] 2139 — 

CII]2326 1.00±0.08 

[NeIV]2425 — 

[Oil] 2472 — 

[MgVII]2510 — 

[MgVII]2629 — 

Hell 2733 — 

[MgV]2783 0.18±0.08 

MgII2800 1.39±0.17 

OIII2836 — 

OIII3133 — 

Hell 3203 — 

[NeV]3346 — 

[NeV]3426 — 

CII] 2326 flux 

(10-i''ergs-icm-2) 3.71±0.30 



Nucleus 



NW EELR 

(1".4 NW) 



4.84±0.20 
1.71+0.11 
0.53±0.12 
0.16+0.04 
0.22+0.04 
0.2.3+0.04 
2.28+0.09 
0.09+0.02 
1.00+0.04 
1.09±0.05 
0.20+0.03 
0.11+0.02 
0.27+0.03 
0.12+0.02 
0.21±0.04 
1.70±0.07 
0.13+0.02 
0.26±0.04 
0.14±0.04 
0.54±0.05 
1.45±0.08 



3.06±0.32 
1.40±0.24 



1.89±0.17 

l.OOiO.08 
0.55+0.09 



0.30+0.14 
1.30±0.17 
0.23±0.17 



<0.67 



7.10±0.25 



4.32±0.34 



Table 6. Emission-line integrated fluxes (corrected for Galactic 
reddening and normalized to the CII] 2326 flux) for the different 
regions in 1303+091. The details of the apertures are in the cap- 
tion of Fig. 10. The quoted errors correspond to the line fitting 
errors. 



observed polarization is on average ^ 8 per cent. In addi- 
tion to this, inspection of the broad-band images of 0850-206 
and 1303-1-091 shows no evidence for a point source in the 
nucleus of cither galaxy. Direct AGN light is therefore not 
an important mechanism in the UV continuum of 0850-206 
or 1303-1-091. Thus, nebular emission, young stars and scat- 
tered AGN light will be the main mechanisms contributing 
to the UV continuum emission of both sources. 

The theoretical nebular continuum for Case B recombi- 
nation (Osterbrock 1989) can be derived from the strength 
of the recombination lines, H/3 in particular. As the ob- 
served wavelength range for 0850-206 and 1303-1-091 did 
not cover the II/3 line, the nebular continuum emission of 
these galaxies was instead derived by using the recombina- 
tion line Hell 1640 and the average line ratio Hell / H/J ~ 
3.18 given by McCarthy (1993) for radio galaxies. The spa- 
tial profiles of the nebular continuum emission for 0850-206 
and 1303-1-091 are shown in Figs. 2 and 7, respectively. To 
calculate these profiles, a slice perpendicular to the disper- 
sion direction, of 40 A width in the observed frame, contain- 
ing the Hell emission line weis extracted from the total two- 
dimensional spectrum of each galaxy. Then, adjacent line- 
free continuum slices were extracted, combined and scaled 
to have the same spectral width as the HcII slice. The com- 
bined continuum slice was then subtracted from the Hell 
slice to produce a continuum-free Hell spatial profile, which 
was used to generate the nebular continuum omission pro- 
file (using the task NEBCONT in DIPSO). For both sources, 
the nebular contribution is generally more important in the 
extended regions than in the nucleus. In 0850-206, essen- 
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Source 


A range (A) 


Po 


(%) 


Pobs {%) 


Fncb (%) 


Fstar (%) 


Fscatt (%) 


0850-206 


1700 


-1850 


6- 


-16 


23±3 


~ 22 


~ 


~ 78 




2900 


-3050 


6- 


-17 


15.6±1.3 


~ 11 


~ 


~ 89 


1303+091 


1700 


-1850 


6- 


-16 


6.2±1.2 


~ 11 


~ 0-50 


~ 39-89 




2900 


-3050 


6- 


-17 


10.1±1.7 


~ 11 


~ 0-29 


~ 60-89 



Table 7. Approximate relative contributions of the nebular continuum stellar population Fgtar and scattered AGN light Fscatt 

to the total continuum emission of 0850—206 and 1303+091 in the rest wavelength ranges 1700 - 1850 A and 2900-3050 A, evaluated for 
a 2-arcsec wide aperture centred at the continuum centroid of the galaxies. Also are given the ranges of intrinsic continuum polarization 

Po for dust scattering, assuming a half-cone opening angle of 45° and angles between the cone axis and the line of sight in the range 
50°— 90° (Manzini & di Serego Alighieri 1996), together with the observed polarization Fobs in the two wavelength ranges. 



tially all of the light south of 2 arcsec from the nucleus is 

nebular omission. For the apertures used in the polarization 
analysis, the contribution of the nebular continuum to the 
total continuum emission in 0850-206 is found to be ~ 22 
per cent for the rest-frame wavelength range 1700-1850 A, 
and ~11 per cent for the rest-frame range 2900-3050 A. 
In the case of 1303+091, the nebular contribution is ~ 11 
per cent for both rest-frame wavelength ranges 1700 - 1850 
A and 2900 - 3050 A. Thus, at least ~ 80 per cent of the UV 
continuum in these galaxies will be a combination of stellar 
light and scattered AGN light ^. It is likely that if the con- 
tinuum polarization of the source is low, the contribution of 
a young stellar population will be important. On the other 
hand, if it is the case that the observed continuum polar- 
ization is high, then the scattered quasar light will probably 
dominate the UV continuum emission of the source. 

The relative contributions of the continuum compo- 
nents can be roughly quantified by assuming that the ob- 
served fractional polarization Poba is diluted by unpolarized 
radiation, which is mainly composed of nebular continuum 
emission Fneb and stellar light Fstar- Thus we can use the 
following equation: 



-t scatt I neb \ ^ star 

where Po is the intrinsic undiluted polarization and Fscatt is 
the scattered AGN flux. Assuming that dust scattering is the 
dominant mechanism for 0850-206 and 1303+091 (see Sec- 
tion 6.1.1 below for a discussion about the scatterers), and 
taking into account that Zubko & Laor (2000) showed that 
polarization is almost independent of optical depth, the cal- 
culations derived by Manzini & di Serego Alighieri (1996) for 
optically thin dust scattering cam be used to obtain a good 
estimate of the intrinsic polarization of the scattered radia- 
tion. Their model uses the geometry indicated by the unified 
schemes for radio sources (Barthel 1989) and Eissumes that 
the nuclear radiation is emitted along two diametrically op- 
posed cones with a half-opening angle of 45° . Depending on 
the angle between the cone axis and the line of sight, the ob- 
served intrinsic polarization for a given wavelength will vary. 
Using equation (1) and assuming that the observer's viewing 

t The stellar light will include both the old stellar population 
of the galaxy and any young stellar population. Based on K- 
magnitudes at these redshifts (K-z relation; e.g. McCarthy 1993), 
the old stellar population is likely to contribute ^ 10 per cent of 
the light at 3000 A, and will be negligible at wavelengths < 2500 A. 



angle relative to the cone axis varies between 50° and 90°, 

the approximate contributions of the different components 
to the contirmum emission of 0850-206 and 1303+091 can 
be estimated for the rest-wavelengths ranges 1700 -1850 A 
and 2900 -3050 A. The results are listed in Table 7. 

In the case of 0850-206 the contribution of young stellar 
population is negligible or null, and the scattered AGN light 
dominates the continuum emission, with a relative contribu- 
tion of ~ 80 - 90 per cent depending on the nebular contribu- 
tion for the different wavelength ranges. Based on the results 
obtained for observer's viewing angles of 50° and 90° using 
equation (1), it is likely that the angle between the radio 
axis and the line of sight for 0850-206 is close to 90° , given 
that for angles below 90° the resulting unpolarized contri- 
bution is negative (particularly in the 1700 - 1850 A range). 
This orientation is also consistent with the radio image of 
this galaxy (see Fig. 1). 

On the other hand, in the case of 1303+091, the stellar 
population can account for up to 50 per cent of the contin- 
uum emission if the radio axis is on the plane of the sky, or 
its contribution could be null if the angle between the line 
of sight and the radio axis is close to 50° . 

An alternative analysis consists of fitting the spectra 
of both galaxies with four components: nebular continuum, 
scattered quasar light, young stellar population and old stel- 
lar population. The quasar spectrum used for the fits was a 
radio-loud quasar composite (from Tclfcr ct al. 2002; with 
a power- law index = -0.67 for rest-frame wavelenghts 
A > 1300 A, and derived using over 200 spectra taken with 
the Hubble Space Telescope, with redshifts z>0.33). The 
quasar spectrum was then scattered using Zubko & Laor 
(2000) dust scattering models, both optically thin and op- 
tically thick (see Section 6.1.1). Both old and young stel- 
lar populations were modelled using the 2000 version of the 
GISSEL codes of Bruzual & Chariot (cf. Bruzual & Chariot 
1993, 2003), assuming a Salpeter IMF, solar metallicity and 
instantaneous burst models of age 4 Gyr and 10 Myr (ap- 
proximate age of the radio source) , respectively. The nebular 
continuum contribution was fixed by our observations and 
the other three components were allowed to vary. Another 
variable was included in the fit to account for any plausible 
intrinsic reddening in each galaxy. 

The result of the fits is shown in Fig. 11. The best fit 
in both cases was obtained when the quasar spectrum had 
been scattered by optically thin dust, and for both galEixies 
a small intrinsic reddening was required: Ejj*'^° = 0.087 for 
0850-206, Eg'*'{? = 0.020 for 1303+091. This intrinsic red- 
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Figure 11. Four-component fit to the spectra of 0850—206 (left) and 1303+091 (right). The components of the fit include; nebular 
continuum (dotted line), dust-scattered radio- loud quasar light (dashed line), young stellar population (dot-dash-dot line) and old stellar 
population (thin solid line). The total fit is indicated by a thick solid line. Both the 0850-206 and 1303-f091 spectra have been smoothed 
into 4 A (rest-frame) bins. 



dening could plausibly be due to the same dust that is caus- 
ing the scattering. The contribution of old stellar population 
in the best-fitting model is consistent at K-band wavelengths 
with typical K-band magnitudes for z~1.4 radio galaxies 
(sec K-z relation; McCarthy 1993). The fits arc good in 
general apart from the wavelength region ^ 2000-2200 A 
where the fit underpredicts the observed continuum emis- 
sion (for further discussion see Section 6.1.1). The results 
from the fitting show that in 0850-206 the scattered quasar 
light dominates the UV contirmum emission and there is no 
requirement at all for a young stellar component; whilst in 
the case of 1303-1-091, both scattered quasar light and young 
stellar population have important contributions to the UV 
light: in the wavelength range 1700 -1850 A contributions 
are comparable, while in the 2900 -3050 A range scattered 
quasar light contributes more than 50 per cent. These re- 
sults are in complete agreement with our calculation above 
(see Table 7). 

Characteristic photospheric absorption features of OB 
stars, such as SV 1502 and NIV 1720, were searched for in 
the spectra of both sources, especially in 1303-1-091, where 
the contribution of a young stellar population is suspected to 
be significant. No clear evidence for any of these absorption 
lines was found in any of the spectra, although the S/N 
of the spectra is insufficient for their absence to rule out a 
young stellar population. 



6.1.1 The nature of the scatterers 

The nature of the scattering material in high-redshift radio 
galaxies is an important question, still under debate. Both 
electrons (Fabian 1989) and dust grains (Tadhunter et al. 
1988) have been proposed as possible scattering agents; but 
dust is generally preferred, mainly because of the greater 
scattering efficiency of dust grains over electrons. Polariza- 
tion observations of radio galaxies show that for the implied 
scattered light to be fully scattered by electrons, the re- 
quired mass of ionized gas is very large (^ 10^'^ - 10^^ Mq; 
e.g. di Serego Alighieri, Cimatti & Fosbury 1994; Cimatti 
et al. 1996). Thus, although some light might be scattered 



by electrons (mainly in the central regions), the scattering 
process is likely to be dominated by the dust grains (pro- 
vided that some dust is present). 

Zubko & Laor (2000) present detailed calculations of 
the scattering and polarization properties of light scattered 
by both optically thin and optically thick dust, using the 
Mathis, Rumpl & Nordsieck (1977) Galactic dust model. 
They find that the spectrum of the scattered light is quite 
different in the two cases. In the optically thin case the emer- 
gent scattered spectrum is bluer than the incident one for 
A > 2200 A, and it has a broad dip centred at A ~ 1500 A. 
In the optically thick case the emergent spectrum is gen- 
erally gray scattered for A > 2200 A, and shows a drop at 
A < 2200 A. Although the total scattered intensity depends 
on optical depth, the polarization of the scattered light 
varies with wavelength in a very similar way for both op- 
tically thin and optically thick cases. The polarization tends 
to decrease slightly from 10000 A to ~ 2000 A, shows a dip 
centred at ~ 1800 A, and then tends to increase sharply to- 
wards the far-UV. 

In the case of quasar light scattered by electrons, the 
emergent spectrum will be similar to the incident one, 
since the Thomson scattering cross-section is independent of 
wavelength. Thus, for A > 2200 A, optically thick dust scat- 
tering can mimic electron scattering. 

In Fig. 12 contirmum polarization properties of 0850- 
206 (left) and 1303-1-091 (right) are presented. In each panel, 
from top to bottom the following arc plotted: the total ffux 
spectrum, the observed fractional polarization, the polar- 
ized flux with a composite radio-loud quasar spectrum (from 
Telfer et al. 2002) superimposed, and the ratio between the 
polarized flux spectrum and the composite quasar spectrum. 
The dust-scattering predictions from Zubko & Laor (2000) 
are overplotted on this bottom plot (see caption for details). 

This composite quasar spectrum can be considered to 
be an average quasar spectrum, and thus it can be compared 
with the polarized ffux spectra of 0850-206 and 1303+091 
(Fig. 12). It is found that, for both gala^xies, the polarized 
spectrum is slightly redder than the incident one, which is 
in agreement with spectropolarimetric studies of other high- 
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Figure 12. Polarization properties of 0850-206 (left) and 1303+091 (right). For each panel, from top to bottom: the unsmoothed total 
flux spectrum in an expanded scale to show the continuum, the observed fractional polarization, the polarized flux with a composite 
radio-loud quasar spectrum [from Telfer et al. (2002), normalized to the polarized flux at rest- frame 2000 A] superimposed, and the 
ratio between the polarized flux spectrum and the composite quasar spectrum with the Zubko & Laor (2000) models overplotted [for 
optically-thin (thin lines) and optically-thick (thick lines) dust cases, with scattering angles of 90° (solid lines) and 120° (dashed lines)] . 
Filled circles and crosses indicate respectively continuum bins and emission line bins with their underlying continuum. 



redshift radio galaxies (e.g. Cimatti et al. 1996; Kishimoto 
et al. 2001). The ratio between polarized and incident spec- 
tra can be directly compared with the scattering efficiency 
X intrinsic polarization [of. equation (1)]; the dust scatter- 
ing predictions of Zubko & Laor (2000) were used for this 
comparison. It can be seen that the broad dip shortward 
of ~ 2200 A (rest-frame) predicted by both optically-thin 
and optically-thick dust models is found in both galaxies, 
more pronounced in 1303-1-091. It therefore appears that 
dust dominates the scattering process in both 0850-206 and 
1303+091. 

The detection of features in the scattered light spec- 
trum around 2200 A is particularly interesting when com- 
pared to the extinction laws of different galaxies. The ex- 
tinction curve of the Milky Way shows a strong 2200 A dust 
feature; this feature is loss pronounced for the Large Magel- 
lanic Cloud (LMC), is even weaker for the Small Magellanic 
Cloud (SMC), and is entirely absent in nearby starburst 
galaxies (e.g. Calzetti, Kinney & Storchi-Bergmann 1994). 
These differences are likely to be related to the metallicity 
of the galaxy and the composition of the dust. The fact that 
a 2200 A dust feature is detected in the total flux spectra 
(Fig. 11) and even more strongly in the polarized flux spec- 
tra (Fig. 12) of the two radio galaxies studied in this paper, 



implies that the dust in these distant galaxies is difl'erent in 
nature or composition to that of nearby starburst galaxies, 
and is much more comparable to Galactic dust. On the other 
hand, the flts to the total continuum spectra (Fig. 11) show 
that, although the observed continuum and the model have 
the same shape around in the wavelength range ~2000- 
2200 A, the model (which uses a Galactic dust scattering 
prediction) slightly underpredicts the continuum emission 
in this range. This suggests that the 2200 A feature is not as 
strong as in the models, and that the dust iu these galaxies, 
although similar, may not be exactly the same as Galactic 
dust. 

6.2 Ionization of the emission-line regions 

In the above sections where the results are described, the 
spatial variation of the line ratio CIII]1909 / CII]2326 along 
the position of the slit for 0850-206 and 1303+091 were 
presented. For 1303+091 it was found that the ionization 
state peaks in the nuclear region, with a possible second 
peak at ^ 1.2 arcscc to the north-west of the nucleus. For 
0850-206 the highest ionization state was found ^ 0.5 arcsec 
offset of the nucleus to the north, then it decreases at both 
sides of this location. 
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Figure 13. Diagnostic diagrams for different regions in 0850— 
206 (triangles) and 1303+091 (stars). Bigger symbols indicate 
the nucleus, smaller ones denote EELR. Solid lines represent op- 
tically thick single-slab power-law (F^oc;/™) photoionization mod- 
els (MAPPINGS) with a=-1.0,-1.5 and a sequence in ioniza- 
tion parameter (2.5xlO~^<U<10~i). The dot-dash-dot line cor- 
responds to photoionization models including matter-bounded 
clouds (10^'^<Apj,i/i<10; Bincttc, Wilson & Storclii-Bcrgmann 
1996). Shock (dashed grid) and shock+prccursor (dotted grid) 
models are included (Dopita & Sutherland 1996); each se- 
quence corresponds to a fixed magnetic parameter (B/y^=0,1,2,4 
//Gcm~^/^) and changing shock velocity (200<i!<500 kms~^). 



In order to study in more detail the ionization state 
of the emission-hne gas in the two galaxies, the line ra- 
tios (derived from the total intensity spectra) for the diflfer- 
ent regions have been compared with theoretical predictions 
of several ionization models. The resulting diagnostic dia- 
grams arc presented in Fig. 13, which involve the following 
line ratios: CUT] 1909 / Hell 1640, GUI] 1909 / CII]2326, and 
[NeIV]2425 / [NeV]3426. The points correspond to the dif- 
ferent emission-line regions defined along the slit (for details 
of the apertures see captions of Figs. 5 and 10). 

The data are compared with the following ionization 
models: (i) optically-think power-law photoionization, with 
spectral indexes a = 1.0 and -1.5 (solid linos); (ii) pho- 
toionization including matter-bounded clouds (dot-dash-dot 
line); (iii) shock- ionization, with shock velocities in the 
range 200<t;<500 kms~^ (dashed grid); and (iv) shock- 
ionization including a photoionized precursor (dotted grid). 
For details about how the different models were generated, 
see Solorzano-Inarrea et al. (2001) and references therein. 

By comparing the three diagnostic diagrams, it can be 
seen that the points of the nuclear regions (bigger symbols) 
of the two galaxies fall close together in the three diagrams, 
indicating that they have a similar state of ionization. Both 
nuclei are well reproduced by the photoionization predic- 
tions with a ~ -1.0 and U ~ 0.025. Note that they could also 
be explained by the mixed-medium photoionization models, 
since these models can be slightly tuned to better fit the 
data. The shock-ionization predictions, however, clearly fail 
to reproduce the ionization state of the nuclear regions in 
the two galaxies. 

The EELR of 1303+091 (open stars) are well ex- 
plained by the photoionization models with a~-1.0, al- 
though the data can also be consistent with low velocity 
shock-ionization predictions. Note, however, that in the mid- 
dle and lower diagrams there is no data for the SE region, 
since none of the neon lines was measurable in this region 
and therefore the ratio [NelV] / [NeV] could not be calcu- 
lated. Only a lower limit could be derived for the same line 
ratio for the NW region. Thus, the ionization of the EELR 
in 1303+091 cannot be well constrained. 

The EELR of 0850-206 (triangles) fall well apart from 
each other in the diagrams. Whilst the ionization of the 
northern region N can be explained by a ~-1.5 photoionizar 
tion predictions and by shock+precursor models with shock 
velocities of ~500kms~^; the southern region SI is best 
explained by mixed-medium photoionization models with 
0.05 < Am/i < 0.1. On the other hand, the ionization of the 
outermost southern region S2 is poorly constrained and thus 
can be reproduced by cither the photoionization models or 
shock-ionization predictions with or without a precursor gas. 

In summary, it is found that the ionization state of the 
nuclear regions in the two galaxies can only be reproduced 
by photoionization predictions. The photoionization models 
give also reasonable fits to the data of the extended regions. 
However, their ionization is less constrained, and with the 
exception of the SI region in 0850-206, whose ionization 
is only consistent with the mixed-medium photoionization 
models, the ionization of the other EELR could also be ex- 
plained by shock-ionization predictions. In particular, the 
SE region of 1303+091, which coincides with the location of 
a radio knot, was found to have the lowest ionization state 
and the highest emission linewidths; these features are ex- 
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pected when the gas is perturbed by strong interactions with 
the radio source structures (e.g. Clark et al. 1998). 



6.3 Variations with radio source size 

Clearly, continuum studies of only two sources are not 
enough to draw significant conclusions about the evolution 
of the continuum properties of radio galaxies with radio 
source age. Nevertheless, the results presented in this pa- 
per are in line with what might have been expected. 

If a jet-induced starburst occurs when the radio 
hotspots pass through the host galaxy, by a few x 10^ years 
later (typical age of a radio source) the starburst luminos- 
ity will have decreased to ~ 10 per cent of its peak value 
(Best, Longair & Rottgering 1996). Therefore, the presence 
of a young stellar population, possibly jet-induced, may di- 
lute the scattered AGN light in smaller (younger) radio 
sources and, consequently, the observed continuum polariza- 
tion would be lower. On the other hand, for larger (older) 
radio sources, in which the contribution of a starburst to 
the total optical flux will be dramatically decreased, the 
scattered AGN light dominates the UV continuum emission. 
Our findings are consistent with this scenario. 0850-206, the 
larger radio source of the two (Drad ~ 118 kpc), presents a 
high continuum polarization, averaging ~ 17 per cent; scat- 
tered quasar light dominates the UV continuum emission 
and the starlight contribution is null or negligible. By con- 
trast, 1303-1-091 (Drad ~ 73 kpc) presents a lower continuum 
polarization, averaging ~ 8 per cent, and in this galaxy the 
young stellar population could account for up to ~ 50 per 
cent of the UV continuum emission. 

These results can be combined with those of previous 
polarization studies of high-redshift radio galaxies, in order 
to explore if there is any correlation between continuum po- 
larization and size of the radio source. Fig. 14 shows the 
fractional continuum polarization versus radio source size 
for a sample of high-redshift radio galaxies taken from the 
hterature (see caption of the figure for details). It must be 
emphasized that this set of radio galaxies consists of several 
different subsamples, many of which are biased (e.g. towards 
high polarization objects). In addition, these samples are at 
different rcdshifts, making the wavelength range studied dif- 
ferent and bringing evolution effects into play, thus making 
the objects not directly comparable. However, even despite 
this, it is intriguing that there is a hint of a weak correla- 
tion between polarization and radio source size (Spearmann- 
Rank test gives a 97 per cent significance level): smaller radio 
sources do not show high polarizations, and the larger the 
radio source the higher its polarization can be. 

The relation between radio source size and polarization 
will be better tested when our full sample of radio galax- 
ies at the same rodsliift (z~1.4), with a large range of ra- 
dio sizes, has been observed. Also, it will be interesting to 
sum the spectra of the smaller radio sources in the full sam- 
ple, with similar characteristics to 1303+091, to search at 
higher S/N for photospheric absorption features of OB stars, 
which would directly prove the presence of a young stellar 
population. 
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Figure 14. Observed fractional continuum polarization of high- 
redshift radio galaxies versus their radio source size. Diflferent 
symbols indicate different redshift ranges: 0.5<z<1.0 (filled 
stars), 1.0<z<2.0 (filled triangles), 2.0<z<3.0 (filled squares) 
and 3.0<z<4.0 (filled circles). Also plotted are the two ra- 
dio galaxies from our sample presented in this paper (big open 
squares). The data corresponding to the other galaxies have been 
taken from di Scrcgo Alighicri ct al. (1989, 1993, 1996); Tadhuntcr 
et al. (1992); Cimatti et al. (1993, 1996, 1997); Dcy & Spinrad 

(1996) ; Dcy ct al. (1997) and Vernet ct al. (2001). Information 
about the sizes of radio sources has been obtained from Dunlop 
et al. (1989); iVIcCarthy, van Breugel & Kapahi (1991); Rottgering 
et al. (1994); CariUi, Owen & Harris (1994); CariUi et al. 

(1997) ; Cohen et al. (1997) and Kapahi et al. (1998). 



7 SUMMARY AND CONCLUSIONS 

This paper presents the results obtained from VLT spec- 
tropolarimetric observations of two powerful radio galaxies 
at z~ 1.4. Analysis of their scattered flux spectra and their 
polarization properties indicates that the scattering process 
in these galaxies is dominated by dust, whose 2200 A fea- 
ture suggests that it is similar to Galactic dust in nature. 
The larger radio source (0850-206) presents a contirmum 
fractional polarization reaching as high as ^ 24 per cent at 
~ 2000 A (rest-frame). For this galaxy, scattered AGN Ught 
must dominate the UV continuum, the nebular emission con- 
tributes up to ~ 22 per cent and there is no requirement for 
any young starlight contribution. The smaller radio source 
(1303-1-091) presents a lower continuum polarization averag- 
ing ~ 8 per cent across the observed wavelength range, and 
its nebular continuum contribution to the UV continuum is 
~ 11 per cent. For this galaxy, the starlight could contribute 
between and 50 per cent, depending on the radio axis ori- 
entation, with multi-component fit to the spectral energy 
distribution favouring the upper end of this range. In both 
galaxies, the position angle of the electric vector is almost 
independent of wavelength and within 15° of perpendicular 
to the radio axis. 

The emission-line properties of the galaxies are also 
analysed. Illumination by the central AGN is likely to be the 
dominant ionization mechanism of the nuclear and EELR in 
both galaxies, although the region coincident with the radio 
knot in 1303-1-091 shows features reminiscent of jet-cloud 
interactions. 
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A compilation of polarimetric observations from the lit- 
erature shows hints of a correlation between continuum po- 
larization and radio source size, albeit for incomplete and 
biased samples. This is broadly in line with what might have 
been expected; that in smaller radio sources the presence of 
a young stellar population (possibly jet-induced) may dilute 
the scattered AGN light, while in larger (older) radio sources 
the presence of young stars is negligible and the scattered 
AGN light dominates the UV continuum. It is clearly neces- 
sary to observe and analyse the rest of our sample to confirm 
this result with a complete and unbiased sample, in order to 
allow significant conclusions to be drawn about the evolu- 
tion of the continuum properties of radio galaxies with radio 
source size. 
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